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The adsorption of fluorine, chlorine, bromine, and iodine diatomic molecules on graphene has been 
investigated using density functional theory with taking into account nonlocal correlation effects by 
means of vdW-DF approach. It is shown that the van der Waals interaction plays a crucial role in the 
formation of chemical bonding between graphene and halogen molecules, and is therefore important 
for a proper description of adsorption in this system. In-plane orientation of the molecules has been 
found to be more stable than the orientation perpendicular to the graphene layer. In the cases of F2, 
Br2 and I2 we also found an ionic contribution to the binding energy, slowly vanishing with distance. 
Analysis of the electronic structure shows that ionic interaction arises due to the charge transfer 
from graphene to the molecules. Furthermore, we found that the increase of impurity concentration 
leads to the conduction band formation in graphene due to interaction between halogen molecules. 
In addition, graphite intercalation by halogen molecules has been investigated. In the presence 
of halogen molecules the binding between graphite layers becomes significantly weaker, which is 
in accordance with the results of recent experiments on sonochemical exfoliation of intercalated 
graphite. 



PACS numbers: 73.22. Pr, 73.20.Hb 

I. INTRODUCTION 

A monolayer of graphite, commonly known as 
graphene, the first truly two-dimensional crystal (onc- 
atom-thick), which became experimentally available in 
the last yearsi Today, graphene is at the focus of many 
research activities worldwide. Remarkable electronic 
properties of graphene make this material a promising 
candidate for a large variety of electronic applications^ 

There have been a number of publications focused on 
theoretical investigation of covalent-, ionic- and metal- 
bonded impurities on graphene (for example, see RefsQ - 
@) ■ These studies are primarily based on the density func- 
tional theory (DFT) combined with standard functionals 
for electronic exchange-correlation effects, such as local 
density approximation (LDA) and generalized gradient 
approximation (GGA). The employment of (semi-)local 
approximations looks quite reasonable so far as the sys- 
tems with either covalent or metallic bonding are con- 
cerned. However, in the case of weakly bounded systems, 
such as molecules or molecular compounds, dispersion 
forces play a significant role. Such interactions are essen- 
tially nonlocal and therefore cannot be properly taken 
into account using local or semi-local density function- 
als. 

In this work, we will focus on molecular impurities 
on graphene. In particular, we consider diatomic halo- 
gen molecules on an ideal graphene lattice. The mo- 
tivation for this choice is the following. First, these 
molecules themselves are covalently bonded and repre- 
sent an example of realistic weakly reactive adsorbates. 



Second, because of the relative chemical inertness of the 
molecules, the interaction between these molecules and 
the graphene sheet is rather weak and cannot be prop- 
erly described within the electronic structure methods 
widely used nowadays in condensed matter physics (such 
as conventional LDA or GGA approximations in the den- 
sity functional theories). Accurate calculations of the ad- 
sorption energies, which are necessary to judge how the 
impurities can be removed from the surface, are beyond 
these standard approaches which make their calculation 
a challenging problem. Last but not least, among the 
fundamental aspects of studying impurities on graphene, 
in some cases molccularly doped graphene represents it- 
self an unavoidable outcome of its production processes. 
Recently Windekvist et al& have reported a method of 
graphene fabrication based on sonochemical exfoliation 
of bromine-intercalated graphite. It was shown that ul- 
trasonic treatment of graphite with absorbed bromine 
molecules results in intensive graphene-fiakes formation. 
It is worth mentioning that intercalated graphite by itself 
has many potential applications^ 

In this paper we show that the adsorption of halogens 
on graphene is determined generally by the van der Waals 
(vdW) interactions. Besides the vdW interaction we also 
found a considerable ionic contribution for F2 and I2 ad- 
sorbates. We demonstrate that the ionic component of 
the interaction arises due to the electronic transfer from 
graphene to the molecules. Furthermore, the ionic inter- 
action gives rise to a long-range tail of adsorption curves, 
which determines the adsorption at large distances. 

In this work we also apply a van der Waals density 



functional approach (vdW-DF) 10 ! 11 to investigate struc- 
tural and energetic properties of graphite intercalated by 
different diatomic halogen molecules. We show that the 
interlayer distance and binding energies in graphite are 
also dependent on the exchange-correlation functional. 
Using vdW-DF method we have calculated interplanar 
binding of graphite layers in the presence of halogen 
molecules, which was found to be significantly weaker 
than the binding in pure graphite. 

The rest of the paper is organized as follows. In Sec. II 
we briefly describe methods of the calculations. In Sec. 
Ill we present the results of adsorption of F2, CI2, Br2, 
and I2 on graphene. Section IV is devoted to the analy- 
sis of electronic structure and ionic interaction between 
graphene and its adsorbatcs. The results of graphite in- 
tercalation by halogen molecules and corresponding dis- 
cussion are given in Sec. V. In Sec. VI we briefly sum- 
marize our results. 
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FIG. 1: Schematic representation of different high-symmetry 
adsorption sites for diatomic molecules on graphene surface 
(top view). In the case of in-plane orientation of the molecules 
the center of mass is located above: (a) a bridge site (the mid- 
dle of C-C bond) along z-axis, (b) a hollow site along s-axis, 
(c) a hollow site along y-axis, (d) a hollow site along ^x-axis, 
(e) a bridge site along y-axis. Molecules perpendicular to the 
surface are located above: (f) a bridge site, (g) a hollow site, 
(h) a top site. 



II. COMPUTATIONAL DETAILS 
A. Crystal structure 

The structure of the systems under consideration was 
modelled using a supercell approach. The supercell con- 
sists of one graphene layer with a molecule on the top and 
a vacuum region to avoid spurious interaction between 
periodic images of the supercell in z-direction. In order 
to examine the effects of impurity concentration two su- 
percells containing 24 and 32 carbon atoms per single 
molecule were considered. The corresponding graphene 
layers were constructed from rectangular (3 x 2V3)a and 
(4 x 2>/3)a unit cells, respectively. The height of the su- 
pcrcclls perpendicular to the surface was chosen to be 40 
A. Investigating the supercells of such a size allows us to 
judge how the interactions between molecules can affect 
the system properties. 

We used the lattice constant of graphene equal to 
a=2.459 A in accordance with experimentally obtained 
value for graphite at low temperatures^ The distortion 
of graphene layer in the presence of adsorbates was not 
taken into account. Since the interaction between the 
molecules and graphene is pretty weak the distortions 
of carbon lattice is supposed to be small, anyway. In- 
teratomic distances in the halogen molecules were also 
taken from the experimental dataJ^ 

Both in-plane and perpendicular-to-plane orientations 
of the molecules on graphene have been investigated. We 
considered five in-plane and three perpendicular-to-plane 
configurations corresponding to high-symmetry adsorp- 
tion sites as illustrated in FigfTJ 

B. DFT computational details 

Ground-state energies and electronic density distribu- 
tions have been calculated using the plane-wave pseu- 



dopotential method as implemented in the Quantum- 
ESPRESSO simulation package ^ 4 ' 15 Exchange and cor- 
relation effects have been generally taken into account 
using vdW-DF nonlocal functional (see the next subsec- 
tion for details). In order to obtain some comparative 
results the GGA functional in Pcdrcw-Burkc-Ernzcrhof 
(PBE) parametrizatior*i& also was used. 

In our calculations we employed an energy cutoff of 30 
Ry for the plane-wave basis and 300 Ry for the charge 
density. Self-consistent calculations of the Kohn-Sham 
equations were carried out employing the convergence 
criterion of 10 -8 Ry. For accurate Brillouin-zone inte- 
gration the tetrahedron scheme^ and (16 x 16 x 1) 
Monkhorst-Pack k-point mesbi^ were used. 

Geometry optimization algorithms were not used in 
our study. Atoms in graphene layer and interatomic dis- 
tances in the molecules were assumed to be fixed. In 
order to find equilibrium distances of adsorbates we per- 
formed a set of explicit energy calculations for different 
separations from the surface. The most stable positions 
of the molecules were determined comparing the energies 
of different high-symmetry adsorption sites. 



C. Non-local corrections 

In the recent years several reasonable methods have 
been developed to improve the performance of DFT 
in description of sparse matter J£r— Most of them are 
based on the adiabatic connection fluctuation-dissipation 
framework (ACFD), which is a computationally very de- 
manding approach and, in practice, is limited to small 
systems only . 22 i 23 

Here, to calculate adsorption energies, we use an ap- 
proach proposed by Dion et ali^ The main idea of the 
method is that the ACFD correlation energy can be sig- 
nificantly simplified by applying a series of reasonable 
approximations resulting in the expression dependent on 
the electronic density only. Although being developed 
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only recently, up to now this method has been widely 
verified on a large variety of examples and shows the 
transferability across a broad spectrum of interactions, 
such as ionic, covalent, and vdW3^ 

In the vdW-DF method, the exchange-correlation en- 
ergy functional consists of several contributions, 

E xc [n] = E^ PBE [n] + E LDA [n] + E?[n), (1) 

where the first term corresponds to the exchange part 
of the revised PBE (revPBE) functional^ E^ DA is the 
LDA correlation energy, and E™ l [n] is the nonlocal cor- 
relation correction, which is calculated in the following 
way: 

Ef = i J d 3 rd 3 r'n(r)(f>(r,r')n(r') 7 (2) 

where n(r) is the electronic density and </>(r, r') is a func- 
tion incorporating many-body density response (for de- 
tails see Reffiol). It should be emphasized that in out- 
work we evaluate the nonlocal correction (Eqj2j in a per- 
turbative way, i.e. using only GGA-based (semi-local) 
electronic density distribution. This choice looks quite 
reasonable since it has been previously shown that the 
effects due to the lack of self-consistency are negligible*^ 

III. ADSORPTION OF MOLECULES ON 
GRAPHENE 

As we have already mentioned, we consider two differ- 
ent sizes of the superccll in order to figure out the influ- 
ence of intermolecular binding on the investigated prop- 
erties. We found that the adsorption energy at higher 
concentrations of halogens is slightly lower, however, this 
difference docs not exceed 10 meV and could be con- 
sidered as insignificant. Due to this fact we do not ex- 
pect any considerable changes in adsorption for different 
molecular concentrations. In the rest of this section we 
present the results only for the smaller supercell contain- 
ing 24 carbon atoms. 

We found the in-plane bridge position (FigQJi) to be 
the most stable configuration for all the molecules under 
consideration. The binding energies for this particular 
geometry are shown in Fig(2] as a function of the dis- 
tance between the surface and the center of mass of the 
molecule. 

It is clear from Fig J2] that the account of the vdW in- 
teraction gives a considerable contribution to the binding 
energy for each molecule examined. One can see that the 
larger the radii of atoms in the molecule the stronger is 
the vdW binding to the surface. This fact can be easily 
figured out using the interpretation of the vdW interac- 
tion as an interaction of induced dipoles. 

One can see from Fig(5] that the binding energy of 
the molecules decreases quite slowly with the distance. 
This fact is the most evident for the fluorine molecule 
(Figl2k), due to the relatively small energy scale. More- 
over, the behavior of fluorine binding with graphene at 



large distances is fairly close to inverse proportionality 
and, therefore, indicates the presence of a strong ionic 
interaction. Indeed, at the distance of several angstroms 
the electronic density of graphene does not overlap with 
that of adsorbed molecules, which means that covalent 
interaction cannot occur in this case. Furthermore, the 
slowly vanishing tail of adsorption curves remains, even 
without taking into account the vdW correction, which 
confirms the ionic nature of such long-range interaction. 

In the case of originally neutral constituents of the sys- 
tem, the ionic interaction could emerge only due to the 
charge transfer between them. Comparison of the bind- 
ing energy for different halogens at relatively large dis- 
tances (more than 6 A) shows that in the absence of the 
vdW interaction the binding of graphene with fluorine 
is stronger than with other molecules. Assuming that 
this interaction is essentially ionic, we expect the larger 
charge transfer in F2 case. By the same reasoning the 
charge imbalance for I2 should be larger than for both 
Br2 and CI2. A more detailed discussion on the ionic 
interaction will be given in the next section. 

For comparison we also show in FigJ5] the results ob- 
tained by the standard semi-local GGA functional in 
PBE parametrization. While the equilibrium distances of 
vdW-DF and PBE approaches agree reasonably well, the 
adsorption energies in the case of PBE are highly overes- 
timated. The latter fact is not surprising since GGA-likc 
functionals do not incorporate nonlocal correlation ef- 
fects, which are responsible for the vdW interaction. It 
should be also pointed out that PBE functional yields 
the wrong asymptote for the adsorption curves with the 
exception of fluorine, where the vdW contribution at the 
large distances is not as large as the ionic one. 

Adsorption energies for different orientations of the 
molecules are summarized in Table HI For both in-plane 
and perpendicular-to-plane orientations the less stable 
configurations correspond to the hollow adsorption site 
(FigfTp,c,d,g). This result is quite reasonable as long as 
the vdW interaction dominates in the system. Indeed, 
at the same distance the strength of the vdW interaction 
is determined by electronic densities corresponding to in- 
teracting components of the system*^ It turns out that 
the hollow site corresponds to the minimum of electronic 
density in graphene, which accounts for the lower vdW 
interaction in this case. 

The difference between various adsorption energies of 
in-plane (or perpendicular-to-plane) orientations of the 
molecules is small. This is especially clear for the ro- 
tated molecules at the hollow adsorption site (FigQ~p,c,d). 
In contrast, the adsorption energies of in-plane orien- 
tations are sufficiently different from the energies of 
perpendicular-to-plane orientations, exceeding thermal 
energy at room temperature (^25 meV). 

We assume that the molecular bond length cannot be 
significantly changed during the adsorption process be- 
cause the obtained molecule-surface interaction is much 
weaker than the intramolecular binding in halogens J£ 
Due to the same reason we do not expect the dissocia- 




FIG. 2: Adsorption energies of halogen molecules on graphene as a function of the separation between graphene and the center 
of mass of the molecules: (a) F2, (b) Cb, (c) Br2, (d) I2. Blue, green and red curves correspond to the total binding energy 
calculated within vdW-DF approach, nonlocal vdW correction (Eq[2]). and to the binding energies excluding vdW correction, 
respectively. Brown curve corresponds to the PBE binding energy. Adsorption curves are shown only for the most stable 
molecular configuration (Fig[TJi). 



TABLE I: Binding energies in graphene-halogen system for different high-symmetry adsorption sites of the molecules, and 
corresponding equilibrium distances (in parentheses). Adsorption sites denoted as in FigfT] Values are given in meV and A, 
respectively. 





in-plane 




perpendicular-to- plane 


a 


bed 


e 


f g h 



F 2 -231 -227 -225 -226 -222 -214 -187 -218 

(3.17) (3.23) (3.24) (3.23) (3.26) (3.47) (3.70) (3.46) 

Cl 2 -259 -248 -249 -250 -251 -208 -193 -210 

(3.58) (3.65) (3.64) (3.65) (3.63) (4.21) (4.34) (4.21) 

Br 2 -296 -289 -287 -285 -288 -229 -207 -214 

(3.74) (3.77) (3.77) (3.77) (3.78) (4.51) (4.67) (4.45) 

I 2 -523 -514 -508 -501 -512 -375 -352 -379 

(3.80) (3.84) (3.83) (3.83) (3.81) (4.74) (4.85) (4.73) 



tive adsorption of halogens on graphene, which usually 
occurs for light molecules on metallic surfaces^ We note, 
however, that a careful description of this phenomenon 
requires atomic dynamics to be taken into account, and 
therefore should be investigated separately. 



IV. ELECTRONIC STRUCTURE AND IONIC 
INTERACTION 

In the previous section we have shown that the vdW 
interaction gives a significant contribution to the ad- 
sorption energy, and, therefore, plays an important role 
in formation of bonding between halogen molecules and 
graphene. However, besides the vdW part there is an- 
other source of interaction which gives rise to the bond- 
ing. In particular, this contribution especially important 
for fluorine and iodine adsorbatcs. In order to clarify the 
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FIG. 3: Total electronic density of states of halogen molecules 
adsorbed on graphene calculated for different supercells. Left 
(right) panel corresponds to the supercell containing 24 (32) 
carbon atoms. Vertical line accentuates the Fermi energy. 
Densities of states are shown for the respective ground state 
distances obtained without vdW corrections (see FigHJ. 



nature of such an interaction let us analyze the density 
of electronic states in the absence of the vdW forces. 

In Fig|3]we have plotted the density of states (DOS) 
for each considered molecule in two different supercells, 
which is calculated for the lowest-energy structural state 
without taking into account nonlocal correction to cor- 
relation energy (excluding third term in EqJT]). In the 
vicinity of the Fermi level DOS exhibits resonances which 
are entirely localized on p-orbitals of halogen molecules. 
The main differences between the spectra for different 
molecular concentrations are related to the width of im- 
purity bands. In both cases the impurity states close to 
the Fermi level remain unoccupied. 

Analysis of the electronic bands shows that the impu- 
rity states close to the Fermi level correspond only to 
one electron per spin, i.e. to one electronic level in non 
spin-polarized case. This means that the observed im- 
purity band at high concentrations of molecules cannot 
be associated with the hybridization between molecular 
and graphene orbitals. In fact, broadening of the impu- 
rity level caused by the hybridization between molecular 



p^-orbitals from different images of the supercell form- 
ing parabolic-like DOS, which is very similar to that of 
a one-dimensional chain. Indeed, for the most energet- 
ically favorable configuration (FigOJi) the distance be- 
tween molecules in x-direction of the supercell is smaller 
than the distance in other directions. As the distance 
between atoms in the molecule increases, the separation 
between different molecules becomes smaller, which leads 
to hybridization of molecular orbitals and broadening of 
the electronic levels. This also accounts for the absence 
of broad impurity bands for lower molecular concentra- 
tions where the lateral interactions between molecules are 
negligible. The shape of the band for the smallest super- 
cell (24 carbon atoms) clearly shows formation of one- 
dimensional band, with typical divergencies (van Hove 
singularities) at the edges. The width of this band, that 
is, the distance between the singularities, is estimated as 
twice the effective hopping parameter. 

For low concentrations, the impurity states are located 
well below or above the conical (Dirac) point, which is 
easily recognizable in the density of states. However, this 
is not the case for Br2 and I2 adsorbed at their highest 
concentrations (FigJ5}:,d). In this situation, the Dirac 
point lies inside the impurity band. Moreover, since the 
impurity states for Br2 and I2 are broadened enough, 
they can be considered as conduction bands leading to a 
metallic behavior of the system. Therefore, the conduc- 
tivity of graphene can be adjusted varying the concen- 
tration of Bi - 2 and I2 adsorbates. There are preliminary 
experimental evidences that doping by iodine can essen- 
tially increase the conductivity of graphene^ 

As can be seen from Fig|3]thc Fermi level is shifted be- 
low the Dirac point for the all molecules except chlorine. 
A shift downward means that the electrons are donated 
by graphene to F2 , Br2 and I2 molecules (acceptors) . As 
discussed above, the electronic transfer is not surprising 
for the investigated systems. It turns out that the as- 
sumption of an ionic interaction between graphene and 
its adsorbates, given in the previous section, is in a good 
agreement with the peculiarities of electronic structure. 

For a more detailed description of the ionic interaction, 
let us analyze the charge transfer phenomenon. As the 
molecules come closer to graphene, redistribution of the 
electronic density results in the formation of an interface 
dipole layer, which is responsible for the ionic interac- 
tion in the system. In FigUwe show the plane-averaged 
density differences An(z) calculated along z-axis of the 
supercell for each investigating molecule at the distance 
of 5 A from the surface. The most noticeable feature 
of FigfJ] is that there are almost no changes in electron 
density for chlorine in comparison with other molecules. 
This seems to be a reasonable result because the Fermi 
energy coincides in this case with the Dirac point. The 
electron redistribution An(z) for the other molecules has 
nearly the same shape, which is fairly close to the shape 
of a typical dipolar distribution. An(z) exhibits a main 
minimum (maximum) corresponding to localization of 
the charge near the carbon layer (molecule). In the in- 
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FIG. 4: In-plane-averaged electron density difference calcu- 
lated along z-axis of the supercell for F2, CI2, Br2, and I2 
molecules on graphene. z — corresponds to position of the 
carbon layer. In the inset an integrated amount of the charge 
transferred from the surface to the molecules as a function of 
the distance is displayed. 

set of FigQ] we also show an integrated charge, trans- 
ferred from graphene to the molecules as a function of 
the distance between them^ As expected, the amount of 
the transferred charge decreases as the distance between 
graphene and the molecules increases. As can be seen 
from FigfJl at the equilibrium distance fluorine molecule 
accepts more electrons than either I2 or Br2. This fact 
completely agrees with our previous results and explains 
the large ionic interaction between fluorine molecule and 
graphene. 



V. INTERCALATION OF GRAPHITE 

It is known that the van der Waals interaction plays 
an important role in understanding and proper descrip- 
tion of the interplanar bonding in graphite^r— In spite 
of the fact that the standard local density approxima- 
tion (LDA) completely misses dispersion interaction (at 
least for long distances), it reproduces correctly the in- 
tcrlayer spacing in graphite and even the binding energy 
in some cases 1^ This surprising result seems, however, 
to be accidental because it is impossible to reproduce 
accurately the experimental compressibility of graphitic 
structures within the LDA<22 Recently it was shown 3 - 3 - 
that the vdW-DF approach yields not only the proper 
distance between graphite sheets, but also the binding 
energy, which is very close to the experimental estima- 
tions available toda y 34 ' 35 The vdW-DF method was also 
successfully applied to the investigation of energetics in 
potassium-intercalated graphite yielding a good agree- 
ment with experiment^ 

In this section we use the vdW-DF approach to de- 
scribe structural and energetic properties of graphite in- 



TABLE II: Equilibrium interlayer distances, absorption ener- 
gies, and interlayer binding energies for the pure and halogen- 
intercalated graphite. Absorption and binding energies are 
given in meV per absorbed molecule and in meV per carbon 
atoms in the supercell, respectively. 





C24 


C24F2 


C24CI2 


C 2 4Br 2 


C24I2 


d,A 


3.6 


6.0 


7.0 


7.5 


7.4 


AE aba , meV/mol 





648 


605 


542 


114 


AE bind , meV/C 


-49.5 


-12.9 


-13.6 


-14.6 


-22.9 



tcrcalated by halogen molecules. We consider graphite in 
j4B-stacking containing 48 carbon atoms and 2 halogen 
molecules per unit cell. Absorbed molecules are assumed 
to be uniformly distributed over graphite occupying in- 
plane bridge sites in accordance with the most stable po- 
sition on graphene layer (FigH^,). Other parameters of 
the calculations were taken to be the same as described 
in Sec.IIB. 

It was recently demonstrated that graphene sheets are 
impermeable even to light gases including heliumi 37 ' 38 
For this reason, we cannot expect that the intercalation of 
halogen molecules in graphite is a result of a penetration 
of the molecules underneath the top surface layer. On the 
contrary, we suppose that the formation of intercalated 
graphite involves a consecutive exfoliation of the whole 
graphite structure. 

In Table [TT] we summarize equilibrium distances, ab- 
sorption energies, and interlayer binding energies cal- 
culated for the investigated systems^ The interlayer 
distance in intercalated graphite is approximately twice 
larger than in the pure case. We note, however, that the 
estimation of the interlayer spacing is not very accurate 
since we ignore in our calculations a distortion of the 
intercalated graphite lattice, which will inevitably take 
place in the real situation. The absorption energies given 
in Table [TT] are positive because their calculation involves 
a positive energy of the exfoliation process. It should be 
pointed out that among the other considered molecules, 
iodine molecule has a minimal absorption energy value. 

The binding energies obtained for the pure graphite 
are in a good agreement with the previous vdW-DF 
calculations! 33 ' 36 Comparing the interlayer binding for 
the pure and for the intercalated graphite, we conclude 
that the binding becomes significantly weaker in the pres- 
ence of halogen molecules. Moreover, the binding energy 
is lower for the more massive molecules because of their 
stronger interaction with the carbon layers, as discussed 
in Sec. III. Our results are in a qualitative agreement with 
the previous LDA calculations for bromine-intercalated 
graphite £ Although the interlayer binding for our par- 
ticular impurity concentration seems to be overestimated 
(in absolute values) in comparison with the LDA results, 
the general tendency is reproduced quite well and al- 
lows to confirm the interpretation of the experiment on 
sonochemical exfoliation of intercalated graphite given in 
R,efi. 
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VI. CONCLUSION 

We have performed a first-principles investigation 
of the adsorption of diatomic halogen molecules on 
graphene. We have shown that a major part of the 
binding energy in the system corresponds to van der 
Waals interaction, and, therefore, cannot be properly 
described by means of standard (semi-)local approxima- 
tions to exchange-correlation energy. It has been found 
that even at the large distances (up to 10 A) there is 
a non-zero interaction between graphene and F2, Br2, 
I2 molecules. This interactions are of ionic nature and 
arise due to the charge transfer from graphene to the 
molecules. 

In contrast to F 2 and Cl 2 adsorbates, the densities 
of states for graphene in the presence of Br 2 and I 2 
molecules exhibit impurity bands right above the Fermi 
level, which can result, for large enough concentration, 
in the metallization of graphene. Moreover, the width of 
the impurity band is directly dependent on the impurity 



concentration. 

Finally, investigating properties of intercalated 
graphite we have found that interlaycr binding becomes 
significantly weaker than in pure graphite. This is 
in agreement with the previous experimental and 
theoretical investigations of intercalated graphite. 
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